On the stability of general relativistic geometric thin disks 
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The stability of general relativistic thin disks is investigated under a general first order pertur- 
bation of the energy momentum tensor. In particular, we consider temporal, radial and azimuthal 
"test matter" perturbations of the quantities involved on the plane z = 0. We study the thin disks 
generated by applying the "displace, cut and reflect" method, usually known as the image method, 
to the Schwarzschild metric in isotropic coordinates and to the Chazy-Curzon metric and the Zipoy- 
Voorhees metric (7-metric) in Weyl coordinates. In the case of the isotropic Schwarzschild thin disk, 
where a radial pressure is present to support the gravitational attraction, the disk is stable and the 
perturbation favors the formation of rings. Also, we found the expected result that the thin disk 
models generated by the Chazy-Curzon and Zipoy-Voorhees metric with only azimuthal pressure 
are not stable under a general first order perturbation. 
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I. INTRODUCTION 

In the last decades, analytical axially symmetric disk 
solutions have appeared in both Newtonian and General 
Relativistic formulations. In the context of General Rel- 
ativity, several exact solutions were found, among them 
the static disks without radial pressure studied by Bon- 
nor and Sackfield Q and Morgan and Morgan 2\. Disks 
with radial pressure and with radial tension have been 
considered by Morgan and Morgan Q and Gonzales and 
Letelier 0, respectively. Also, stationary disk models 
including electric fields [5|, magnetic fields [6|, and both 
electric and magnetic fields have been studied. Self 
similar static disks were considered by Lynden-Bell and 
Pineault [8| and Lemos |9(. Furthermore, the superpo- 
sition of black holes with static disks were analyzed by 
Lemos and Letelier ITTl IT^ and Klein Bicak 
and Ledvinka [Llj found relativistic counter-rotating thin 
disks as sources of Kerr type metrics, and Bicak, Lynden- 
Bell and Katz ^| obtained static disks as sources of 
known vacuum spacetimes from the Chazy-Curzon met- 
ric [lgL IT?! and Zipoy-Voorhees metric 0, fTsj [ . also 
Bicak, Lynden-Bell and Pichon [2fJ found an infinite 
number of new static solutions. Recently, exact solutions 
for thin disks with single and composite halos of mat- 
ter [2l]], thin disks made of charged dust thin disks 
made of charged perfect fluid [23, and thick disks [24[ 
were obtained. For a survey on relativistic gravitating 
disks, see |25j . 

It is well known that the formation of stellar systems 
lays in its stability. The study of stability, analytically 
or numerically, is vital to the acceptance and applica- 
bility of the different models found in the literature to 
describe stellar structures found in Nature. On the other 
hand, the study of different types of perturbations, when 
applied to stellar structures, usually give an insight on 
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the formation of bars, rings or different stellar patterns. 
Also, a perturbation may cause the collapse of an stable 
object with the posterior appearance of a different kind 
of structure. 

In Newtonian theory, perturbations have been made 
extensively through the years in analytical calculations 
and numerical experiments. An analytical treatment of 
the stability of thin disks can be found in Refs. [2(| [27j 
and references therein. In General Relativity, the sta- 
bility analysis is usually done studying the particle mo- 
tion along geodesies and not perturbing the energy mo- 
mentum tensor of the fluid and its conservation equa- 
tions. The stability of particle motion along geodesies 
have been studied transforming the Rayleigh criteria of 
stability [2^, [2|j into a General Relativistic formulation, 
see |30| and reference therein. By using this criterion, the 
stability of orbits around black holes surrounded by disks, 
rings and multipolar fields were analyzed in [3(j . Also, 
it was employed in fH) to test stability in the isotropic 
Schwarzschild thin disk, and thin disks of single and com- 
posite halos generated from the Buchdahl solution [3l| 
and Narlikar, Patwardhan and Vaidya solution [3^ by 
using the displace, cut and reflect method. For the sta- 
bility of circular orbits in stationary axisymmetric space- 
times see Ref. [33l ESI • The stability of circular orbits of 
the Lemos-Letelier solution f° r the superposition of a 
black hole and a flat ring, as well as other as pec ts of this 
solution, are considered in detail in |35l l36t |3jj . Bicak, 
Lynden-Bell and Katz 0] analyzed the stability of the 
Chazy-Curzon and Zipoy-Voorhees thin disks without ra- 
dial pressure or tension studying their velocity curves and 
specific angular momentum, finding that they are not sta- 
ble for highly relativistic disks. A general stability study 
of a relativistic fluid, with both bulk and dynamical vis- 
cosity, perturbing its energy momentum tensor was done 
by Seguin [38j. However, he consider the coefficients of 
the perturbed variables as constants, i.e., local perturba- 
tions. Fact that in general is too restrictive. 

The purpose of this work is to study numerically the 
stability of several analytical thin disk solutions in the 



2 



performing a general first order perturbation in the tem- 
poral, radial and azimuthal components of the quanti- 
ties involved in the energy momentum tensor of the fluid 
and analyzing the corresponding perturbed conservation 
equations of motion. The perturbations considered do 
not modified the background metric obtained from the 
solution of Einstein equations, i.e. they are treated as 
"test matter" . 

The article is organized as follows. In Section II-A 
we present the perturbed conservation equations for the 
general case, when the energy momentum tensor is di- 
agonal with the energy density and the three pressures 
or tensions all different from zero. Later, in Section II- 
B, we specialize these perturbed equations for the thin 
disk case where we have a distributed energy momen- 
tum tensor with support on the plane 2 = 0. In Sec- 
tion III, we analyze the perturbed conservation equations 
to study the stability of the isotropic Schwarzschild thin 
disk, which contains radial and azimuthal pressures. In 
Section IV, we analyze the perturbed conservation equa- 
tions to study the stability of the Chazy-Curzon and 
Zipoy-Voorhees thin disks with only azimuthal pressure. 
It is worth to mention that different Chazy-Curzon and 
Zipoy-Voorhees thin disks models with radial tensions or 
pressures can be built by using a different set of coor- 
dinates, but these models do not usually satisfy all the 
physical requirements, see 0- Finally, in Section V, we 
summarize and discuss the presented results. 



where U^, X^, and are the 4- vectors of the or- 
thonormal tetrad 

tf" = e-* 1 (l,0,0,0), 
X* = e -*3 (0,1,0,0), 

y = — (0,0,1,0), 

r 

z n = e -*3( ; o,0,l), (3) 

which satisfy, 

-u^u„ = x^x^ = y% = z"z M = l, 

U»X M = EP% = U<*Z„ = 

= X^Z„ = Y^Z^ = 0. (4) 

The quantity is the time-like four- velocity of the fluid 
and the quantities X^, Y^ and Z^ are the space-like 
principal directions of the fluid. 

In the axially symmetric case, all the quantities in- 
volved in the energy momentum tensor are functions of 
the coordinates (r, z) only, and so are the coefficients of 
the perturbed conservation equations. With this in mind, 
we construct the general perturbation A p of a quantity 
in the form 



II. PERTURBED EQUATIONS 

A. General case 

Let us consider the general static, axially symmetric 
metric 



A p (t, r, tp, z) = A»(r, z) + 5A»(t, r, <p, z), 



(5) 



where A^ (r, z) is the unperturbed quantity and 
<5A M (i, r, tp, z) is the perturbation. Applying the pertur- 
bation (JSJ in J2J) we obtain that 



ds 2 = -e^'dt 2 + e 2 ^r 2 d^ 2 + e 2 ^(dr 2 + dz 2 ), (1) 

where Vfi, \&2 and ^3 are functions of the variables (r, z). 
(Our conventions are: G = c = 1. Metric signature 
+2. Partial and covariant derivatives with respect to the 
coordinate x M denoted by , p, and respectively.) 

The energy momentum tensor of the fluid sat- 
isfies Einstein field equations G^ = kT^ v , and, when 
written in its rest frame, is diagonal with compo- 
nents {—p,p r ,P<p,Pz), where p is the energy density and 
{pr,P<p,Pz) are the radial, azimuthal and axial pressures 
or tensions respectively. With these definitions, the en- 
ergy momentum tensor can be written as 

= pU^U v + p r X^X u + p v Y^Y v + p z Z»Z' y , (2) 
I 



Tp" (t,r,tp,z) = T>» (r,z) + 5T»» (t,r,tp,z). (6) 

Henceforth, we assume that the perturbed energy mo- 
mentum tensor does not modify the background metric 
found by solving the Einstein field equations G M „ = kT^, 
i.e. the perturbation ST^ acts as a test fluid. With this 
assumption, the perturbed energy momentum equations 
for the problem can be written as 



{8T^), y = 0. (7) 

Writing in explicit form the four equations from JjJ and 
discarding terms of order greater or equal to (S) 2 , we ob- 
tain the perturbed conservation equations for the system 



/! = t 
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+SU r F(t, r, pU 1 ) + SU z F(t, z, pU 1 ) + OT*F(i, r, Pr X r ) + SZ*F(t, z,p z Z z ) = 0, (8) 

p = r 

SU^pU 1 ) + SX: r ( Pr X r ) + SX%( Pr X r ) + 5Y^( Pv Y*) + 5Z^ z { Pz Z z ) + Sp rir (X r X r ) + 28U t (pU t Tl t ) 
+26X r G{r, r, Pr X r ) + 5X z F{r, z, Pr X r ) + 26Y* (p^T^) + 5Z r F(r, z lPz Z z ) + 25Z z (p z Z z F r zz ) 
+Sp(U t U t T r u ) + 6 Pr G(r, r, X r X r ) + 8p v {Y^T r ¥V ) + 6 Pz (Z z Z z F r zz ) = 0, (9) 

fjb = ip 

SU^pU 1 ) + 5Xf r (p r X r ) + 5Y*(p v Y*) + 5Y a a { Pv Y*) + 5Z%(p z Z z ) + 5p v , v (Y*Y<?) 

+SX v F(tp, r,p r X r ) + SY r F&, r, Pip Y^) + SY Z F{ V , z,p v Y*) + 5Z*F(<p, z, Pz Z z ) = 0, (10) 

p = z 

SU%{pU*) + SX z r ( Pr X r ) + 5Y z v ( Pv Y<?) + SZ z z ( Pz Z z ) + SZ%( Pz Z z ) + 5p z , z {Z*Z z ) + 25U 1 {pU^) 
+2SX r ( Pr X r F z rr ) + 5X z F(z,r,p r X r ) + 28Y* '(P^T^) + 5Z r F{z, r,p z Z z ) + 2SZ z G(z, z, Pz Z z ) 
+S P (U t U t F z t ) + S Pr (X r X r F z r ) + S Pv (Y^F z vv ) + S Pz G(z, z, Z Z Z Z ) = 0, (11) 

where 



F{I,J,K)=K,j + K(2r I IJ + T%j), 
G(I,J,K)=K„j + K(T I IJ 



(12) 
(13) 



and Tp are the Christoffel symbols. 

We want the perturbed tetrad to remain orthonormal. 
Therefore, to guarantee the orthonormal conditions Q 
we obtain the following relations, 





= 5X r = SY V : 


= 5Z Z 




8X l 


= HiSU r , 


Zi = 


-Xr/Ut, 


8Y l : 


= 




-YJU U 


5Z* - 


= &SU", 


e 3 = 


-z z /u u 


5X V 


= U5Y\ 


u = 


—X r /Yp, 


6X Z 


= &sz r , 


& = 


—X r JZ z , 


5Y Z 


^HeSZ^, 


£ 6 = 


-Y v /Z z . 



It is clear from Eqs. (|14fl that in order to have a con- 
sistent perturbation model, the tetrad ought to be per- 
turbed. In our case, we are interested in perturbations in 
the velocity components and for that reason the t com- 
ponent of the tetrad must be perturbed. 

B. Thin disk case 

Thin disks are characterized by a distribution valued 
energy momentum tensor that is a delta function with 
support on the plane z = 0. In this case, the axial pres- 
sure p z is equal to zero and the energy density, radial and 
azimuthal pressures or tensions are only functions of the 
radial coordinate r. So, the energy momentum tensor 



T" v = {<jU^U v + p r X"X v + Pip Y^Y u )8{z). (15) 

where a is the surface energy density and (U^, X^, Y^) 
are the components of the tetrad previously defined. 
Note that when we define the energy momentum tensor 
as above, the fluid quantities (cr,p r ,p v ) are not the "true" 
fluid quantities {o , ,p r ,p~ tp ). These two sets of quantities 
are related by 

& = Vfh^cr, p~ r = y/g^p r , p~ v = yfgTzP^- (16) 

Using the perturbed conservation equations and 
the distributed energy momentum tensor l|15fl we find 
that 



{5T^),J(z) + ST" v [6{z)l v = 0. (17) 

In this work, the thin disk metrics considered are found 
by appl ying th e well known displace, cut and reflect 
method (15L |21| . For that reason, the metric components 
depend on the radial coordinate r and the absolute value 
of z. If we define the derivative of the absolute value 
of z in z = equal to zero (|z|. z =o = 0), and perform 
in (I17fl the integration on the coordinate z (J ^/g zz dz), 
we obtain that the second term is equal to zero and the 
perturbed conservation equations for thin disks can be 
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(5T^). iU \ z = = 0. (18) 

In general, we are interested in perturbations of the 
four-velocity and the thermodynamic quantities on the 
disk. Thus, we set the perturbations on the Z M tetrad 
vector equal to zero. Also, the perturbation 5X V is not 
related to the perturbation in the radial or azimuthal 
velocities and we can set it equal to zero. Note that 
perturbations in the tetrad are independent of the per- 
turbations in the fluid quantities and the assumptions 
made above do not change the equation of state of the 
fluid. So, we have 

SZ f = 5Z r = SZ V = §X V = 0, (19) 
which implies from (|14|l that 



8U Z = SX Z = SY Z = SY r = 0, (20) 

and the perturbed tetrad is still orthonormal. 

Due to the form of the perturbed energy momentum 
equation for thin disks l|15(l . the connection coefficients 
are only functions of the radial coordinate. Therefore, all 
the coefficients on the perturbed conservation equations 
depend only on r and we construct a general perturbation 
of the form 



SA"(t, r, tp) = 5A^(r)e l{k ^ wt) . (21) 

Hereafter 5 A* = 5A^(r). With the perturbation JgTJl 
and conditions l|14(l and l|19f) . the perturbed conservation 
equations (|18|l can be written in explicit form as 



/! = t 



\i = r 



8ir r (aU* + ZiP r X r ) + <5C/ r [F(i,r, ( j[/ t ) + ^ rPr X r + ^F{t,r, Pr X r )} 
+<5^[ifc v (CT[/ t +6jvr p )] + Sai-iwU'U 1 ) = 0, 



5p r , r (X r X r ) + SU^-iwiaU* + £ip r X r )] + 5cr(l7*l7*r&) 
+8p r G{r,r,X r X r )+5 Pv {Y^T r w ) = 0, 



(22) 

(23) 
(24) 



where we set Sp z — in order to maintain the state 
equation of the fluid on the disk. The above equations 
are the basic equations for thin disks stability and will be 
studied in the next sections. The general equations 
(|llf> are presented for completeness and for future works 
on the stability of thick disks and other static axially 
symmetric structures. 



III. THIN DISKS WITH RADIAL AND 
AZIMUTHAL PRESSURES 



a) , and applying the displace, cut and reflect method to 
the resulting isotropic Schwarzschild metric. With this 
operation, we find a thin disk with surface energy density 
and equal pressures given by |21j . 



16mai?Q 
" 7r(m + 2,R ) 5 ' 

Pr =P V = 



8m 2 aRl 



n(m + 2R ) 5 (m-2R ) 



(26) 
(27) 



We start analyzing the isotropic Schwarzschild thin 
disk. This disk is found by setting the metric functions 
(*i,* 2 ,*3) in as 



= In 



* 2 = *a 



2R-m 



2R + m 



In 



2R 



(25) 



where Ro means to evaluate R at z = 0. Recall that 
we are not using the true fluid quantities, th e q uantities 
lPE|l and lj27|l are related to the ones cited in |21j by Eqs. 
(|16fl . From here and through the rest of the manuscript, 
all the calculations are made in the variables (<y, Pr ,p v ) 
and the final results presented in the figures expressed in 
terms of the true physical variables. 

We want our perturbations to be in accordance with 
the equation of state of the fluid, i.e. p = p(r) and a — 
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FIG. 1: Profiles of the energy density perturbation amplitudes of the isotropic Schwarzschild thin disk with parameters 
(a=0.5,m=0.4). The first three different w modes for the first three wave number k are plotted, we see that increasing 
the w mode the number of oscillations within the disk increases while increasing the wave number k the amplitude decreases. 



dp = Sp r = 8p v = p^ r dr, (28) 
5a = o>dr, (29) 

from which we find the useful relation 

S p = f^j 8a. (30) 

Substituting SU r and SU V in Eq. from Eqs. 
and (J52J, and using relation (|30p> . we find a second order 
differential equation for the energy density perturbation 

c _ _ r _ e 



AsSa irr + B s 8a :r + C s 5a = 0, (31) 

where (As, B$,Cs) are functions of (r, a, m, w, k v ), see 
Appendix A. 

Note that the exact thin disk metrics we are consid- 
ering are infinite in the radial direction. So, in order to 
study the stability of the disks we need a criterion to 
make a cut-off in the radial coordinate to create a finite 
disk. From [2lJ we see that the energy density of the 
disks decreases rapidly enough in principle to define a 
cut-off radius. The cut-off radius R cu t of the disk is set 
by the following criterion: the matter within the thin disk 
formed by the cut-off radius is more than 90% of the total 
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FIG. 2: Profiles of the true pressure perturbation, physical radial velocity and physical azimuthal velocity amplitudes of 
the isotropic Schwarzschild thin disk with parameters (a=0.5,m=0.4). We see that the qualitative aspects of the pressure 
perturbation are the same of the energy density perturbation. The amplitude of the radial velocity increases when we get near 
the edge of the disk, in that circumstances we must compared these velocity values to the particle escape velocity to make the 
perturbation consistent with the model. The azimuthal velocity amplitude is almost constant far from the center of the disk. 



energy density is distributed along the plane z = from 
outside the disk to infinity. This allow us, if necessary, 
to treat the outer 10% of the energy density disk as a 
perturbation in the outmost boundary condition. More- 
over, the values of the parameters (o, m) of the finite 
disk we create must satisfy the dominant energy condi- 
tion. For the infinite disk, this condition imposes that 
5" + p r + Pip > f° r au r - This inequality is valid in all 
the disk if m < a, the lower the value of a for a given 
m the faster the surface energy density decreases. Fur- 
thermore, this condition also implies that no tachyonic 
matter is present in the disk and that p r and p v are in 
fact pressures and not tensions. 



The second order equation (|3 1|) is solved with two 
boundary conditions, one in r ~ and the other in the 
cut-off radius set by the criterion established above. At 
rwOwe set the perturbation to be ss 10% of the unper- 
turbed energy density value and at the edge of the disk 
8tj\r=R cut = 0. The last condition is imposed because 
we want the perturbation to vanish when r tends to the 
outer radius. 

Now, we consider the isotropic Schwarzschild thin disk 
with parameters (a = 0.5, m = 0.4). With these parame- 
ters, the outer radius of the disk is set to r = 4 (approx- 
imately 90% of the energy density is inside the disk) . In 

TTi- _ ITI _i _ l * j i _ m _ _ _ r - j 
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density perturbation 5a for different modes of the pertur- 
bation J2J. We see from Fig. [I] that the energy density 
perturbation profiles are stable and have an oscillatory 
character. When we increase the parameter w the num- 
ber of oscillations within the disk increases and when we 
increase the wave number k the amplitude of the oscilla- 
tions decreases. Note that the amplitudes of the modes 
decay quickly when we increase the value of the wave 
number. From Fig. 2] we see that there is a factor of 
approximately 1CP 3 between the modes (w = l,k = 0) 
and (w = 1, k = 2). 

In Fig. we present the amplitude profiles of the 
true pressure perturbation and the physical radial veloc- 
ity perturbation (SU r = SU r ^g rr ) for the first three w 
modes with k = 0, and the amplitude profile of the phys- 
ical azimuthal velocity perturbation {8U V = SU^ \fg^) 
for the first three w modes with k = 1. The modes of 
the azimuthal velocity perturbation with k = are equal 
to zero. We see from Fig. [5] that the pressure pertur- 
bation has the same qualitative aspects of the density 
perturbation. The azimuthal velocity perturbation am- 
plitude shows an oscillatory behavior but the difference 
between it and the radial velocity perturbation is that 
after the first oscillations the maximum value of the am- 
plitude remains almost constant for the rest of the disk. 
Note that the amplitude of the radial velocity pertur- 
bation increases with the radial coordinate. To discuss 
about this fact, we have to take into account that we are 
using in our numerical calculations a finite disk instead of 
the infinite disk exact solution. In the infinite disk there 
is no place for stars to escape from it but in the finite 
disk this could happen. For our disk to be consistent, 
the values of the radial velocity perturbation within the 
disk can not be greater than the particle escape velocity. 
In first approximation, the escape velocity of our disk is 
calculated using Newtonian mechanics. In the Newto- 
nian limit of General Relativity, we have the well known 
relatio n goo = ?/oo + 2$ and the particle escape velocity 
V e = y2|$| can be written as 

For (a = 0.5, m = 0.4, r — 4) the escape velocity is w 
0.42 and we see from the radial velocity perturbation 
profile of Fig. |3 that the particles within the disk do 
not escape. In the case when the outer radius is equal 
to 5 (approximately 92% of the energy density) and the 
same parameters and perturbations are set, the escape 
velocity is « 0.38 and some modes of the radial velocity 
perturbation violates the consistency of the disk and can 
not exist. In other words, these modes describe a disk of 
a certain cut-off radius with no particles inside. 

The case studied in this section corresponds to an ex- 
treme disk because the radial velocity of the particles at 
the cut-off radius are near the particle escape velocity. 
This example was chosen by didactic reasons to show the 



disk model, but in most of the cases these problems do 
not appear. Recall that the perturbation made to the en- 
ergy density was 10% of its unperturbed initial density. If 
we limited our perturbation to lower values, the radial ve- 
locity perturbation through the disk is below the escape 
velocity. Furthermore, we can define the cut-off radius of 
the disk when the energy density within the disk is lower 
than 90%. With lower values for the center perturba- 
tion and percentage of energy density within the disk, all 
modes tested were stable. The phenomenon of particles 
escaping from the disk is more suitable to occur in the 
mode with (k = 0,w = 1). We see from Fig. Q]and Fig. [21 
that increasing the number of the wave vector k the am- 
plitude of the variables considered decreases. All these 
considerations suggest that the isotropic Schwarzschild 
thin disk is stable under a general first order perturba- 
tion of the form (|21|l and favors the formation of rings. 

When we performed the radial cut in the infinite disk 
in order to create the finite disk, we laid outside a fraction 
of the total energy density distributed up to infinity. We 
can model the presence of this exterior energy density, if 
necessary, by perturbing the outer boundary condition, 
i.e., setting in the outer radius the condition Sa\ r =R cut = 
e where e <C Sa\ r =o- The numerical experiments realized 
show that this outer perturbation does not modify the 
qualitative aspects of our solutions. However, depending 
on the value of the perturbation, the quantitative aspects 
do change, but the disk remains stable. 



IV. THIN DISKS WITH ONLY AZIMUTHAL 
PRESSURE 

A. Chazy-Curzon thin disk 

The Chazy-Curzon thin disk in Weyl coordinates 0] 
has only azimuthal pressure, the absence of a radial pres- 
sure turns this solution rather unphysical. Even though, 
one can argue that the stability of the system can be 
achieved if we have two circular streams of particles mov- 
ing in opposite directions, i.e. the counter-rotating hy- 
pothesis. Furthermore, the dragging of inertial frame of 
a rotating disk does not generate a Weyl type metric 
but the counter-rotating hypothesis solve this problem. 
Recently, evidence of disks made of streams of rotating 
and counter- rotating matter has been found j3^. Also, a 
detailed study of the counter-rotating model for generic 
finite axially symmetric thin disk without radial pressure 
can be found in Ref . ^(J • 

The Chazy-Curzon metric is obtained from with 



m r 
A = 7-, 

2R 4 ' 
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FIG. 3: Profiles of the A/B function present in the solution 
of the energy density perturbation for the thin Chazy-Curzon 
disk. The profiles show a singularity (instability) that de- 
pends on the cut parameter a. For a — (1; 1.4; 1.8) and m = 
1, the singularities are approximately at (1.414; 1.980; 2.546) 
respectively. 



FIG. 4: Profiles of the A/B function present in the solution of 
the energy density perturbation for the thin Zipoy-Voorhees 
disk. As the Chazy-Curzon case, the profiles show a singu- 
larity (instability) that depends on the cut parameters (a, b). 
For b = (2; 3; 4), a = 1 and m = 1, the singularities are 
approximately at (2.193; 2.523; 2.978) respectively. 



*3 = A - $. 



(34) 



The surface energy density and azimuthal pressure of the 
the Chazy-Curzon thin disk, after applying the image 
method to the metric with the above definitions, are 
given by [l5|. 



ma 


1 - 


TO / 


2ttRI 


R^ \ 


m?a 





a 2 


2ttR% 


Rl 



Rl 



a 2(* -A ) 



(35) 



(36) 



where Rq, $0 and Ao are the respective functions evalu- 
ated at z = 0. Since the Chazy-Curzon disk has no radial 
pressure, we have p r = Sp r = in the perturbed equa- 
tions. Substituting SU r and SU^ from Eqs. J23J) and 
(|24|l into Eq. ill'L'l) and using a relation similar to i|30|) 
involving the fluid variables Ij35(l and Ij36|l . we obtain a 
first order differential equation in r for the perturbation 
5a of the form 



In the integral argument, the values of (a,m) must sat- 
isfy the dominant energy condition a + p~ v > for all r if 
we want a fluid made of matter with the usual gravita- 
tional attractive property. This inequality is satisfied if 
m/a < PS 1.30. In the case of the Chazy-Curzon 

disk, the dominant energy condition implies that the ve- 
locity of the counter-rotating streams does not exceed 
the velocity of light. We show in Fig. [3] the profile of the 
integral argument A/B present in the solution l|38(l for 
different value of the parameter a. The singular points 
do not depend on the values of the parameters (777, w, kip) 
but on the value of the cut parameter a. The numerical 
experiments realized for the true surface energy density 
6a show instabilities at the same values of the singulari- 
ties present in the function A/B. As far as we can test, 
the Chazy-Curzon disk does not allow any stable mode 
under perturbations of the form l|21() . Stability can be 
achieved if we let p r — and 6p r ^ 0, but in this case is 
being changed the state equation of the fluid on the disk. 



B. Zipoy-Voorhees thin disk 



ASa , 



B5a = 0, 



where (A, B) are functions of (r, a, 777, w, k^), 
pendix B. The above equation has for solution 



(37) 
see Ap- 



!(A/B\dr 



In this section we study the stability of the Zipoy- 
Voorhees thin disks. These disks are obtained from the 
Zipoy-Voorhees metric, also known as the 7-metric. This 
metric represents, in Weyl coordinates, a uniform rod of 
length b — a at a distance a below the z — plane along 
the negative z axis. The Zipoy-Voorhees metric is ob- 
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$ = 
A = 



In 



b — a 
2m 2 
Jb^af- 



R a + \z\ + a 



■In 



R b + \z\ +b 

(Ra + Rb? - (b - a) 21 



4R a R 



(39) 



where i? 2 = r 2 + (\z\ + a) 2 and R 2 = r 2 + (\z\ + b) 2 , and 



same values of the singularities present in the function 
A/B. Moreover, under perturbation of the form (|21|l . the 
Zipoy-Voorhees disk does not have any stable mode. A 
different treatment of the stability of the Chazy-Curzon 
and Zipoy-Voorhees thin disks |15| . made by analyzing 
their velocity curves and specific angular momentum, 
shows that they are not stable for highly relativistic disks. 



* 3 = A - $. 



(40) 



When (b — a) — m, the above solution leads to the 
Schwarzschild metric and when (b—a) — ► 0, to the Chazy- 
Curzon solution. The thin disk formed by applying the 
displace, cut and reflect method has the surface energy 
density and the azimuthal pressure given by |l5j . 



1 



1 



2ir(b- a) \R a o Rbo 
b 



(41) 



Pip 









m f 




i 










b 


-a V 


m 




{t- 


(b 




a) 2 





O 2(* -Ao) 



(42) 



RbO RaO 



O 2(*o-A ) 



where b > a and R a o, Rbo, &o and Ao are the respective 
quantities evaluated at z = 0. Like the Chazy-Curzon 
metric, these disks do not have radial pressure and, in or- 
der to generate a Weyl type metric, the counter-rotating 
hypothesis is applied. 

Following the same steps as in the Chazy-Curzon met- 
ric, we arrive to a first order differential equation that 
yields a solution of the form In this case, the values 
of the parameters (a, b, m) must satisfy the dominant en- 
ergy condition a + p~ v > 0. To achieve this condition we 
must have Ha. 



1 

< - 
a 2 V bRao 



RbO RaQ 



when 



= b 2 ; 



2/3 



aRbo 



a\2/3 
b 



(43) 



(44) 



In Fig. 0] we plot the profiles of the integral argument 
A/B for the Zipoy-Voorhees disk for different values of 
the parameter b. As the Chazy-Curzon disks, the sin- 
gular points that appear in the profiles do not depend 
on the values of (m, w, k v ) but on the value of the cut 
parameters (a, b). Also, numerical experiments show 



V. CONCLUSIONS 

We study the stability of several thin disk models in 
the context of General Relativity using a general first 
order perturbation. The stability analysis performed in 
this work is better than the stability analysis of particle 
motion along geodesies because we take into account the 
collective behavior of the particles. However, this analy- 
sis can be said to be incomplete because the energy mo- 
mentum perturbation is treated like a test fluid and does 
not alter the background metric. The complete analy- 
sis of the thin disks stability needs to take into account 
the perturbation of the metric. This is a second degree 
of approximation to the stability problem in which the 
emission of gravitational radiation is considered. 

We show that the Chazy-Curzon and Zipoy-Voorhees 
thin disks are not stable under perturbations of the form 
(|21|l . This is due to the fact that there is no radial 
pressure to support the gravitational attraction of the 
disk. This appears to be in contradiction with the ex- 
perimental evidence of stellar systems made of streams 
of rotating and counter-rotating matter |39| . but the 
counter-rotating hypothesis assumes the existence of two 
pressureless streams of matter in circular opposite orbits 
which do not interact, see appendix of Ref. |30j , and this 
is not the case. 

In the case of the isotropic Schwarzschild thin disk, 
where we have radial and azimuthal pressures, we formed 
a finite disk by making a cut-off radius and allowing a 
percentage of the unperturbed energy density within the 
disk. The finiteness of our new disk allows the particles 
to escape from it, so we have to compare the particles 
escape velocity with the velocity perturbation profiles if 
we want to have a self-consistent finite thin disk model 
or if we want to discard some perturbation modes. If 
we lower the value of the initial perturbation and/or the 
percentage of the total energy density inside the disk, all 
the modes are stable. The fluid variables, in the isotropic 
Schwarzschild thin disk, present an oscillatory character 
with the amplitudes vanishing when r approaches the 
outmost radius. In the case of the azimuthal perturba- 
tion, the amplitude is almost constant within the disk. In 
general, when we increase the parameter w the number of 
oscillations increases inside the disk and the amplitudes 
decrease. When we increase the wave number k the val- 
ues of the amplitudes decrease abruptly. We note that 
the perturbation l|21l) made in the isotropic Schwarzschild 
thin disk favors the formation of rings. As expected, the 
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bility of thin disks. 
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by Bi and the coefficients of Eq. (23J| by C l7 e.g., 
the first term in (22) has the coefficient A\ multiplied by 
the factor SU r r , the second term has the coefficient A 2 
multiplied by the factor SU r , etc. 

The explicit form of the above equations are obtained 
substituting the fluid variables (p, a) of the isotropic 
Schwarzschild thin disk. 



APPENDIX A: COEFFICIENTS OF THE 
ISOTROPIC SCHWARZSCHILD THIN DISK 
SECOND ORDER PERTURBATION EQUATION 



APPENDIX B: ARGUMENT (A/B) OF THE 
INTEGRAL (l38l) 



The general form of the functions (As, B$, Cs) appear- 
ing in the second order equation (|31|) are given by 



The general form of the argument (A/B) of the integral 
Eq. (55j l are given by 



A s = Aiai, 

B s = Ai(a 1 . r + a 2 ) + A 2 a x , 

Cs = A 1 a 2 , r + A 2 a 2 + ^3^3 + ^4, 



where otx, a 2 and 013 are 



Oil 



a 2 



a 3 



f, r B t + f(B 4 + B b 



Bo 



B 2 



fBi 
B 2 ' 
fC 2 



(45) 



(46) 



and / = (p t r/cr. r ). In Eqs. (|45|l and i|46|) . we denote 
the coefficients of Eq. (22J by Ai , the coefficients of Eq. 



A = a^a^r + F(t, r, oU*)o>! 

-k^oU 1 + ^1^)02 + wU*U*, 
B = 0-17*0!, (47) 

where the functions u\ and 012 are 



Oil = 



a 2 



U*U*TX t +YVYtTr ( Pv>r /o, r ) 



woXJ 1 



(48) 

w(oW + bptpY*) ■ (49) 
The explicit form of the above equations are obtained 
substituting the particular fluid variables [p v ,o) of the 
Chazy-Curzon or Zipoy-Voorhees thin disks. 
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